Strong spin-orbit interaction and helical hole states in Ge/Si nanowires 
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We study theoretically the low-energy hole states of Ge/Si core/shell nanowires. The low-energy valence 
band is quasidegenerate, formed by two doublets of different orbital angular momenta, and can be controlled via 
the relative shell thickness and via external fields. We find that direct (dipolar) coupling to a moderate electric 
field leads to an unusually large spin-orbit interaction of Rashba type on the order of meV which gives rise to 
pronounced helical states enabling electrical spin control. The system allows for quantum dots and spin qubits 
with energy levels that can vary from nearly zero to several meV, depending on the relative shell thickness. 
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I. INTRODUCTION 

Semiconducting nanowires are subject to intense exper- 
imental effort as promising candidates for single-photon 
sourcesji field-effect transistors^ and programmable circuits i^ 
Progress is being made with both group-IV materials^"— and 
III-V compounds, particularly InAs, where single-electron 
quantum dots^ (QDs) and universal spin-qubit control^ have 
been implemented. Proximity-induced superconductivity was 
demonstrated in these systemSj^ii^ forming a platform for Ma- 
jorana fermionsjii^ 

The nanowires are operated in both the electron^™ (conduc- 
tion band, CB) and hola^ ~^-'°''^ (valence band, VB) regimes. 
While these regimes are similar in the charge sector, holes 
can have many advantages in the spin sector. Due to strong 
spin-orbit interaction (SOI) on an atomic level, the electron 
spin is replaced by an effective spin J = 3/2, and even in 
systems that are inversion symmetric, the spin and momen- 
tum are strongly coupled, enabling efficient hole spin ma- 
nipulation by purely electrical means. Holes, moreover, are 
very sensitive to confinement, which strongly prolongs their 
spin lifetimeSfi^^ Also, VBs possess only one valley at the 
r point, in contrast to the CBs of Ge and Si, which is partic- 
ularly useful for spintronics devices such as spin filters^'* and 
spin qubits.-^ Most recently, spin-selective hole tunneling in 
SiGe nanocrystals was achievedf2^ 

In this paper, we analyze the hole spectrum of Ge/Si 
core/shell nanowires, which combine several useful features. 
The holes are subject to strong confinement in two dimen- 
sions and can be confined down to zero dimension (OD) in 
QDs4i21a^ Ge and Si can be grown nuclear-spin-free, and 
mean free paths around 0.5 jum have been reportedi^ Dur- 
ing growth, the core diameter (~5-100 nm) and shell thick- 
ness (~1-10 nm) can be controlled individually. The VB off- 
set at the interface is large, ~0.5 eV, so that holes accumu- 
late naturally in the corei^i^^ Lack of dopants underpins the 
high mobilities- and the charge coherence seen in proximity- 
induced superconductivity tifi 

We find that the low-energy spectrum in Ge/Si nanowires is 
quasidegenerate, in contrast to typical CBs. Static strain, ad- 
justable via the relative shell thickness, allows lifting of this 
quasidegeneracy, providing a high degree of control. We also 
calculate the spectrum in longitudinal QDs, where this feature 



remains pronounced, which is essential for spin-qubit imple- 
mentation. The nanowires are sensitive to external magnetic 
fields, with g factors that depend on both the field orienta- 
tion and the hole momentum. In particular, we find an addi- 
tional SOI of Rashba type (referred to as direct Rashba SOI, 
DRSOI), which results from a direct dipolar coupling to an ex- 
ternal electric field. This term arises in first order of the multi- 
band perturbation theory, and thus is 10-100 times larger than 
the known Rashba SOI (RSOI) for holes which is a third-order 
effecti^ Moreover, the DRSOI scales linearly in the core di- 
ameter R (while the RSOI is proportional to /? '), so that spin- 
orbit interaction remains strong even in large nanowires. Sim- 
ilarly to the conventional Rashba SOI^ii^ilL^ the DRSOI 
induces helical ground states, but with much larger spin-orbit 
energies (meV range) than in other known semiconductors. 

The paper is organized as follows. In Sec.|II]we introduce 
the unperturbed Hamiltonian for holes inside the Ge core and 
provide its exact, numerical solution. The system is very well 
described by an effective ID Hamiltonian, which we derive 
in Sec. |III1 In Sec. |IV] we include the static strain and find 
a strong dependence of the nanowire spectrum on the rela- 
tive shell thickness. The spectrum of Ge/Si-nanowire-based 
QDs is discussed subsequently (Sec.[V]i. In the main section. 
Sec. I VII we analyze the hole coupling to electric fields and 
compare the DRSOI to the standard RSOI. In this context, we 
also show that Ge/Si nanowires present an outstanding plat- 
form for helical hole states and Majorana fermions. Magnetic 
field effects are discussed in Sec. IVIII followed by our sum- 
mary and final remarks. Sec. I Villi Technical details and addi- 
tional information are appended. 



II. MODEL HAMILTONIAN AND NUMERICAL 
SOLUTION 

In cubic semiconductors, the VB states are well described 
by the Luttinger-Kohn (LK) Hamiltonian^^i^ 
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where /j y ; (in units of K) are the three components of the ef- 
fective electron spin in the VB, m is the bare electron mass, hk 
is the momentum operator, and y\ and js = Q-Ji + 3y3)/5 are 
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FIG. 1. Schematic drawing of the systems studied in this paper. Top: 
Excerpt of a Ge/Si nanowire with core radius R and shell thickness 
Rs - R, where the z axis corresponds to the axis along the wire. The 
nanowires are typically several micrometers in length and can there- 
fore be considered infinitely extended, hosting a ID hole gas inside 
their cores. The surrounding Si shell influences the hole spectrum 
through static strain. Bottom: Quantum dots of (effective) length L 
form when the holes are subject to additional confinement in the z 
direction. This can be realized via gates (Refs. 0, |23, and \2w or, 
in principle, by surrounding the Ge with layers of barrier material 
during growth (Ref.l44l). 



the Luttinger parameters in spherical approximation, which is 
well applicable for Ge (yi - 13.35, js = 5.1 l)i^ In studying 
nanowires [Fig. [T] (top)], the LK Hamiltonian must be sup- 
plemented with the confinement in the transverse directions 
{x-y plane), perpendicular to the wire axis z- Since we are 
interested in the low-energy states, we can add two more sim- 
plifications at this stage. First, since the low-energy states 
are located near the core center, we can assume a potential 
with cylindrical symmetry even though the real system is not 
perfectly symmetric. Second, due to the large VB offset, the 
confinement can be treated as a hard wall. 



V(r) 



0, r<R, 
oo, r > R, 



(2) 



with R as the core radius. Given this confinement, the total 
Hamiltonian //lk+ V commutes with the operator F. = L.+J~, 
where L- = -id^ is the orbital angular momentum along the 
wire axis, so that F~ is a good quantum number and the states 
can be classified accordingly,^--- The system is also time- 
reversal symmetric (Kramers doublets), and due to cylindrical 
symmetry one obtains the same spectrum for the same |F.|. 
This is valid for any circular confinement and does not require 
the assumption of a hard wall. We note that, again in clear 
contrast to the CB case, L. is not conserved in the VB. 

The Hamiltonian separates into 4x4 blocks corresponding 
to given F~. By solving Hlk + V numerically, using an ansatz 
analogous to those in Refs. [3a and |37i we find that the low- 
energy spectrum in the Ge core is formed by two quasidegen- 
erate bands, with F^ - +1 /2 each, where the ground (excited) 
states at k- - are of L, a; (|L;| = 1) type. These four 
(in total) bands are well separated from higher bands, and the 
quasidegeneracy indicates that one can project the problem 
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FIG. 2. Low-energy hole spectrum of a Ge nanowire as a function of 
the longitudinal wave number k-. In the unstrained case, 7 = 0, the 
plot is independent oiR, with ?r l{mR-) ^ 0.76 meV fovR = 10 nm. 
Due to time-reversal invariance and cylindrical symmetry, each line 
is a twofold degeneracy, where red (blue) indicates quantum numbers 
F- = ±1/2 (F- = ±3/2). At fc; = the spectrum is quasidegenerate, 
with the lowest states having L- a (ground states) and |L-| = I 
(excited states) character. Dashed red lines result from the effective 
ID model for the lowest subspace, where k- is treated perturbatively. 
The top figure is a plot of the low-energy sector of a strained system, 
y = 40%, illustrating strong dependence on the Si shell thickness. 



onto this subspace. A plot of the spectrum is shown in Fig. |2] 
(bottom). 



III. EFFECTIVE ID HAMILTONIAN 

The present analysis does not, however, allow us to de- 
rive an effective ID Hamiltonian describing the lowest-energy 
states. For this, we integrate out the transverse motion and 
treat k~ in perturbation theory {k^^R < 1). The four eigen- 
states g^: and e±, coiTesponding to ground and excited states 
for F, = ±1/2 at A:- = 0, serve as the basis states in the effec- 
tive ID Hamiltonian. The subscript refers to the sign of the 
contained spin state |±3/2), since the system at k~ - can be 
separated into two 2x2 spin blocks;-^ details of the calcula- 
tion are described in Appendix IB] Knowledge of g± and e-t, 
with eigenenergies Eg = Q and Eg = A, allows us to include 
the fcj-dependent terms of the LK Hamiltonian. The diagonal 
matrix elements take on the form (g±| //lk \g±) - frk^/(2mg), 
{e±\ //lk k±) = trk\l(2mi,) + A, and the nonzero off-diagonal 
terms are of type {e±\Hui,\g^) - iCk-, with C as a real- 
valued coupling constantJS Summarized in matrix notation, 
this yields 



j-f' 



efif 
LK 



: A+ + A-T^ + Cfc.Ty<Tf, 



(3) 



where A± s h^kfim^ ^ +m/)/4- + A/2, and r,-, cti are the Pauli 
matrices acting on {g,e}, {+,-} (see also Appendix lAb. For 
Ge, the values are A = 0.73 lf/(mR-), C = 7.26 Ti^/imR), 



Mg ^ m/(yi + 2ys) — 0.043 m, and m,, — m/{yi + y,) — 
0.054 m. The eigenspectrum 



7 = 



7 = 0.4 



Eg^,{k-) = A+ + yJAl + C-kl 



kl 



(4) 



nicely reproduces all the key features of the exact solution 
and is added to Fig. |2] for comparison, with good agreement 
foihR< 1. 



IV. STATIC STRAIN 

To the above model one needs to add the effects of static 
strain, since the Si shell (radius Rs) tends to compress the 
Ge lattice. A detailed derivation of the strain field in Ge/Si 
core/shell nanowires will be provided elsewhere; here we 
just quote the results needed to calculate the hole spectrum. 
Coupling is described by the Bir-Pikus Hamiltonian //bp, 
Eq. (ICll l, which for Ge (the spherical approximation applies) 
is of the same form as Eq. ^, with kjkj replaced by the strain 
tensor elements sijM Assuming a stress-free wire surface and 
continuous displacement and stress at the interface, symmetry 
considerations and Newton's second law require e„ - e,,, and 
e„. = ej- = ey- = within the core, so that only terms pro- 
portional to j} contribute. Hence, F. remains a good quantum 
number, [//lk + V + //bp, F,] = 0, which allows us to solve 
the system exactly even in the presence of strain, following the 
same steps as described in Sec. |ll] It is important that these 
exact spectra show that the low-energy states [Fig.|2](bottom)] 
separate even further from the higher bands when the Ge core 
is strained by a Si shell, so that the low-energy sector remains 
energetically well isolated and projection onto this subspace 
is always valid. 

In the ID model, strain leads to a simple rescaling of the 
energy splitting A ^ A + (5(y), where < 6{y) < 30 meV for 
< y < oo, with y = {Rs-R)/R as the relative shell thickness. 
Hence, 6 is independent of the core radius, while A oc R^^. 
We note that A ^ 0.6 meV for a wire of /? = 10 nm, which 
makes this energy scale very small. Therefore the splitting 
can be changed not only via R, but also via Rs- In fact, the 
system can be varied from the quasidegenerate to an electron- 
like regime [Fig.|2](top)], where the L^ - and \L~\ ^ 1 states 
are parabolas. 



V. QUANTUM DOT SPECTRUM 

We analyze this feature in more detail by calculating the 
eigenenergies of Ge/Si-nanowire-based QDs [Fig. [T] (bot- 
tom)]. All steps of this calculation are carefully explained in 
Appendix Id] Remarkably, the variability with Rs also trans- 
fers to the QD levels. Figure [3] shows the spectrum as a func- 
tion of confinement length for a wire with both thin and thick 
shells and plots the energy splitting of the lowest Kramers 
doublets as a function of y. For a negligible shell, the states 
lie so close in energy that additional degeneracies may even 
be observed. With increasing Rs, the QD spectrum changes 
monotonically from the quasidegenerate regime to gaps of 
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FIG. 3. Top: Hole energy spectrum in a nanowire-based QD (Ge/Si 
core/shell, R = 5 nm), for both a thin and a thick shell, as a function 
of confinement length L. Each line corresponds to a Kramers pair, 
and dashed lines represent A for comparison. Bottom: Level splitting 
of the two lowest Ki'amers doublets as a function of relative shell 
thickness y and for different lengths L. Static strain, induced via 
the shell, allows continuous tuning of the energy gap over several 
meV, an attractive feature for spin-qubit applications. For details, 
see Appendix IdI 



several meV, which should, in particular, be useful for im- 
plementing spin qubits. 



VI. DIRECT RASHBA SOI AND HELICAL HOLE STATES 

An electric field E^ applied along x couples directly to the 
charge of the hole via the dipole term 



//ed = -e/i.vJC, 



(5) 



with X - r cos((p) as the carrier position in field direction. For 
holes in the Ge core we expect this energy gradient to have 
sizable effects compared to electron systems, since the low- 
energy band is made of quasidegenerate states of different L^ 
character Moreover, E^ will also couple directly to the spins 
due to the SOI in the VB. Projection of //gd onto the subspace 
yields the effective SOI Hamiltonian 



//i 



DR 



//!" = eE,UT,CT„ 



(6) 



referred to as direct Rashba SOI (DRSOI), characterized by 
the coupling constant U - {g+\ {—x) \e+). The form of Eq. ^ 
still resembles that in the CB case, where dipolar coupling 
cannot modify the spins. However, the additional k^TyCTx term 
in H^ makes the key difference to the CB and accounts for 
the SOI featured in the LK Hamiltonian. Indeed, by diagonal- 
izing H^^ + HoR we find that the DRSOI lifts the twofold de- 
generacy, as plotted in Fig.|4] Surprisingly, the effects closely 
resemble a standard RSOI for holes in a transverse electric 




FIG. 4. Dispersion relation for holes in a Ge nanowire of R = 10 nm, 
negligible shell, and an applied electric field E^ along x, calculated 
from /f=^ + Hdr, Eqs. (O and lU, with Hur as the DRSOI Hamilto- 
nian. Hole bands of lower (higher) energy are plotted blue (red). The 
RSOI is about 100 times smaller than the DRSOI and thus negligi- 
ble. Note that the DRSOI shows qualitatively similar features to the 
standard Rashba SOI with dispersion curves shifted along k- against 
each other. 



field (see discussion below). [Again, this is not the case for 
the CB, where //ed does not lift the degeneracy since spin and 
orbit are decoupled (in leading order).] 

As a consequence, when analyzing the eigenstates of H^ + 
Hdr for their spin properties, we find that an electric field gen- 
erates helical ground states, i.e., holes of opposite spin move 
in opposite directions. Figure|5](top) shows the splitting of the 
lowest band when E^ - 6 V/fim is applied to a typical Ge/Si 
nanowire of 5 nm core radius and 1 .5 nm shell thickness. Even 
though RSOI is absent, the result resembles the typical CB 
spectra considered in previous studies, where Rashba SOI for 
electrons leads to two horizontally shifted parabolas in the E-k 
diagramfi^"i^iii2i2i Moreover, the analogy also holds for the 
spins, which are twisted toward the y direction, perpendicular 
to both the propagation axis z and the field direction x. As 
Fig. H] (bottom) illustrates, (Jy) in the ground state is an anti- 
symmetric function of ^., the characteristic feature of a helical 
mode. We note that {J^} = (7-) - throughout, so that the 
spins are indeed oppositely oriented. The values of |7, | around 
the band minima are >l/2, while the spin-orbit (SO) energy, 
i.e., the difference between band minimum and degeneracy at 
k, = 0, is £50 > 1-0 meV. This value exceeds the reported 
100 //eV in InAs nanowires by a factor of 10 (see also Ap- 
pendix |E|j^i^ and further optimization is definitely possible 
via both the gate voltage and the shell thickness. 

We can understand the qualitative similarity of the DRSOI, 
Eq. ©, and RSOI^ 



Hso = aE^ikyJ: - kjy). 



(7) 



by projecting the latter onto the low-energy subspace, which 
yields 



Hr - Hi 



eff 
SO 



aE^S TxCTz 



(8) 



fox k,R < 1, with 5 - {g+\kyJz\e+). Further information 
on Hso, Hr, and the Rashba coefficient a can be found in 
Appendix 10 This formal analogy of //dr and Hr, Eqs. 
(|6]l and ([8]), immediately implies that Ge/Si nanowires pro- 
vide a promising platform for novel quantum effects based 
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FIG. 5. Top: Splitting of the lowest valence band when an electric 
field £j = 6 V//jm is applied to a Ge/Si nanowire of ^ = 5 nm and 
Rs = 6.5 nm. Ground (excited) hole states are plotted blue (red), 
iiso > 1.0 meV, a large value compared to that for InAs (Refs. @ 
and|40|). and the degeneracy at k- = may be lifted via a magnetic 
field (see Fig. [6](. The conventional RSOI for holes is negligible. 
Bottom: Plot of {/,,) for the above system, where (J^) and (J.) are 
zero throughout. In the ground state, the nanowire carries opposite 
spins in opposite directions with \{Jy)\ > 1/2. 



on Rashba-type SOI 
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A particular advantage of 



the DRSOI, as compared to conventional Rashba SOI, is its 
unusually large strength. While the Rashba term for holes 
arises in third order of multiband perturbation theory and 
thus scales with l/(band gap)-, the DRSOI is a first-order ef- 
fect and therefore much stronger,^" Explicit values for Ge are 
U = 0.15 R, S = 0.36/R, and a ^ -OAnm^e, so that, in 
typical nanowires with R = 5-10 nm, //dr dominates //r by 
one to two orders of magnitude (Appendix |Fl). Moreover, siz- 
able RSOI would require unusually small confinement, since 
//r cx R\ In stark contrast, for DRSOI we find //dr oc R, 
which allows one to realize the desired coupling strengths in 
larger wires as well. The upscaling, however, is limited by the 
associated decrease of level splitting (oc R^^) and of the term 
Ck-Ty(T,, (oc /;-') in Eq. Q. 



VII. MAGNETIC FIELD EFFECTS 

The Kramers degeneracy can be lifted by an external mag- 
netic field B, which couples to the holes in two ways, first, 
via the orbital motion, through the substitution hk — > -//iV + 
eA{r), with A{r) as the vector potential, and second, via the 
Zeeman coupling //^ - IkjubB ■ J, where k is a material 
parameter. For B along z (x), parallel (perpendicular) to the 
wire, the ID Hamiltonian is of the form 

Hb.z = I^bB, (Zi (T, + Z2T,(T, + Z3fc,x,cr,,) , (9) 

Hb,.v = fiaB, [Xia, + X2T,(t, + X^k.Ty) , (10) 

where the real-valued constants Z, (X,) are listed in Eq. dGll l 
of Appendix |G] The results agree with recent experiments. 
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FIG. 6. Top: Hole spectrum of Fig. [5] (top) in the presence of 
fij: = IT. The magnetic field opens a gap of 0.30 meV at k. = 0, 
corresponding to a g factor above 5. Bottom: Plot of the ground 
state spin, {/,) and (Jy), where (J-) = throughout. At energies 
within the gap, the Ge/Si nanowire features helical hole states with 
Eso > 1.0 meV, \k-\ == 90 ^m"', and KJ,,>| > 1/2. 



where the g factors in Ge/Si-nanowire-based QDs (multihole 
regime) were found to vary dramatically with both the orien- 
tation of B and also the QD confinement.i^ia^ In the absence 
of electric fields, the ground state g factor g\\ik,) for B^ along 
the wire turns out to be small for fc; = 0, |^||(0)| - 0.1, and 
increases as |A:-| increases. In contrast, the g factor g±ik~) for 
a perpendicular field B^ is large at k- - 0, |gj.(0)| ^ 6, and de- 
creases as l^-l increases, until g±{k.) eventually changes sign 
at \k^\ R! Q.5/R. We note that these results for the ground 
state cannot be directly compared to experimental results in 
the multihole regime, as the g factors in the excited state al- 
ready show a clearly different dependence on fc. In the pres- 
ence of an electric field E^, the effective g\\ and gj_ at fc; = 
may, to some extent, be tuned by the strength of £.1. 

Detailed analysis of the low-energy Hamiltonian yields the 
result that the combination of magnetic and electric fields al- 
lows for optimal tuning of the energy spectrum. For instance, 
Bx - IT opens a gap of 0.30 meV at k~ = in Fig. |5] (top), 
keeping the spin properties for k~ ^ unaffected. This cor- 
responds to |§j.(0)| - 5.2 and is illustrated in Fig. |6l With 
the Fermi level within the induced gap, the spectrum of Fig. |6] 
presents a promising basis for applications using helical hole 
states. Remarkably, an all-perpendicular setup with, e.g., B^ 
along X and £,. along y, //dr.v - -eEylJTy, leads to an asym- 
metric spectrum where only states with one particular direc- 
tion of motion may be occupied, which moreover provide a 
well-polarized spin along the magnetic field axis. As before, 
this does not require standard RSOI. 



VIII. DISCUSSION 

The low-energy properties found in this work make Ge/Si 
core/shell nanowires promising candidates for applications. 



The dipole-induced formation of helical modes proves useful 
for several reasons. First, the strength and orientation of ex- 
ternally applied electric fields are well controllable via gates. 
Second, the DRSOI scales linearly in R, instead of R\ and 
thicker wires remain operational. Third, the system is sensi- 
tive to magnetic fields, and undesired degeneracies at L = 
may easily be lifted, with |^j_(0)| > 5. Finally, helical modes 
with large ^so and wave numbers kf are achievable using 
moderate electric fields of order V/jum. In Fig. |6] with the 
Fermi level inside the gap opened by the magnetic field, these 
are £50 > 1.0 meV and kp ^ 90 //m^', with \{Jy}\ > 1/2, 
and optimization via both the gate voltage and the Si shell is 
possible. For R - 10 nm and thin shells, due to the quaside- 
generacy at y ^ 0, even small electric fields of ~0. 1 V///m 
are sufficient to form helical states with £50 Z 0.3 meV. 
We note that a strong SOI, tuned via electric fields, was re- 
cently reported for Ge/Si nanowires based on magnetotrans- 
port measurements4i 

The nanowire spectrum can be changed from the quaside- 
generate to an electronlike regime, depending on the shell 
thickness. This moreover holds for QD spectra, so that, 
given the strong response to electric and magnetic fields, 
Ge/Si wires also seem attractive for applications in quan- 
tum information processing, particularly via electric-dipole- 
induced spin resonancei^^ila Finally, when combined with a 
superconductor^^ the DRSOI in these wires provides a useful 
platform for Majorana fermions»iir— 
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Appendix A: Representation of spin matrices 

All results presented in this paper are based on the follow- 
ing representation of the spin-3/2 matrices: 
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The Pauli operators r,- (referring to {g, e}) and cr, (acting on 
{+,-)) are defined as 
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and analogously for cr,. 



Appendix B: Basis states for the effective ID Hamiltonian 

In this appendix we outline the calculation of the basis 
states {g+,g-,e+,e-}. For k~ = 0, each of the 4x4 blocks 
for given quantum number F~ and energy E reduces to two 
2x2 blocks, labeled by + according to the sign of the con- 
tained spin state 1+3/2). In the absence of confinement, using 
an ansatz analogous to those in Refs. |3g and [37l the eigen- 
states to be considered are 



Appendix C: Bir-Pikus Hamiltonian 



Referring to holes, the Bir-Pikus Hamiltonian reads 
Hbp = - « + 4^ X ^" "^ ^ Zl ^"-^'^ 

^ ^ i i 

+ ^{e„.[j,,Jy] + c.p.), (CI) 

where a, b, and d are the deformation potentials, e^ = e,, 
are the strain tensor elements, {A,B] = (AB + BA)/2, and 
"c.p." stands for cyclic permutations i^ For Ge, the deforma- 
tion potentials are b ^ -2.5 eV and d ^ -5.0 eVr^ so that the 
spherical approximation d - y3b applies. The hydrostatic 
deformation potential a accounts for the constant energy shift 
of the VB in the presence of hydrostatic strain, and therefore 
does not contribute to (5(y), i.e., the rescaling of the energy 
gap A. 



„i(F-T3/2)«!. I 



- ^JF,±i/2(hbr)e"^^-'^"^^^ 1 + 1/2) , (Bl) 
= V3/F,T3/2(^ihr)e''<^--^'^'>^ 1+3/2) 
+/K±i/2(^ihr)e'<^=^'/2'^ 1+1/2) , (B2) 






where the 7„(x) are Bessel functions of the first kind, and 



_ 1 

thh.lh = 7 



2mE 



h \^yi+ 2y, 
When confinement is present, the eigenstates read 



<l>f=(r, 



^±''^hh±('"''^) + ^±"'/'^'±('"' 



(B3) 



(B4) 



F F 

where the coefficients a_^\b_^ and the energies E are to be 
found from the boundary condition Oj_X/?, (/>) = 0, resulting in 
the determinant equations 



- JF,T3/2{khhR)JF,±\/2(hhR) 

+3JF-±]/2(khhR)JF,+3/2{klhR)- 



(B5) 



By solving the above equations, we find that for + the lowest 
eigenenergy corresponds to F~ = +1/2, and the second lowest 

+ 1/2 

one to F- - +1/2. The associated eigenstates g± = <!)+ and 
e± = <!)+ for the transverse motion are found by calculat- 
ing the coefficients a^ , b^ , a'^ , and fe^ , respectively, 
and serve as the basis states in the effective ID Hamiltonian. 
Normalization requires 



\g±} 



- I drr I ( 

Jo Jo 



#l^±l =1, 



(B6) 



and analogously for e+ . It turns out that the excited states are 
purely heavy-hole-like, b^' - 0, and we choose the com- 
plex phases such that all coefficients are real, with a^' < 0, 
/7±'^^ > 0, and a;;;'^^ > 0. 



Appendix D: Quantum dot spectrum 

When the quantum dot length L is much larger than the core 
radius R, Fig.[Tl the spectrum can be well approximated using 
the eff'ective Hamiltonian for extended states. In the absence 
of external fields, F.^ remains a good quantum number and the 
Hamiltonian 



t/eff _ 
n,v — 



( n^kj 

iCk, 





-iCK 










n'kl 

2«(g 

iCk, 






—iCk, 

2± + ^ + ^(r) ) 



(Dl) 
here explicitly written out in the basis {g+,e_,g_,e+) for il- 
lustration purposes, is 2 x 2 block diagonal with degenerate 
eigenstates. The subspace {§+,e_) corresponds to F^ - -1/2, 
while {g-, e+] corresponds to F- - +1/2. Aiming at the quan- 
tum dot spectrum, we introduce two complex functions g„{z} 
and e„{z), for which we require 



--i -iCk 



( n-k: 
2in^ 

iCk, g + A + 5(r) 



gn(z) 

e„(z) 



- p 

— J-'ti 



gn(z) 
Sniz) 



(D2) 



The associated set of coupled diff'erential equations reads 

J.2 



= -- — g',[{z) - Ce',jiz) - E„gn(z), 
Zm„ 



(D3) 



= -:^<(z) + Cg:,(z) + [A + 6{j) - E„] e„(z), (D4) 

and in addition we demand g„(0) = e„(0) = gniQ - e„{L) = 
due to hard wall confinement at z = and z - L. When the 
dilTerential equations have been solved, these boundary con- 
ditions finally lead to a determinant equation for the eigenen- 
ergies £„, which can be analyzed numerically. The results are 
plotted in Fig.[3j 



Appendix E: Spin-orbit energy in InAs nanowires 

For electrons in an electric field E-^ along x, the Hamiltonian 
for Rashba SOI is of the form 






aEj^ik^CTy — kyCrS), 



(El) 



where a is the Rashba coefficient in the conduction band (Fg) 
and cr, are the Pauli matrices for spin IflM In the following, 
we use the notation a^ = aE^ for illustration purposes. As- 
suming a nanowire in which the electron moves freely along 
the z direction with efifective mass m*, the Hamiltonian of the 
system becomes 



2 7.2 



//'= 



n^k: 



Ini 



— + a^k-CTy, 



(E2) 



with eigenspectrum 



2m* 

n 



S2 



m a^ 



= ^ (^c ± rs'o) - Esc 



(E3) 



The spin-orbit length is defined as /so = fi^/ (m* \a^\), and the 
SO energy, the energy difference between the band minima 
and the degeneracy at k~ - 0, is ^so - m*a\l{2tP'), so that 



^' 2 



(E4) 



We can use Eq. (IE4l i to calculate the spin-orbit energy for InAs 
wires, where /so has recently been measuredi^i^ Using /so - 
127 nm and m* ^ '«buik ~ 0.023 mr- the SO energy in InAs 
is ^so - lOOjueV. Further experiments confirmed that /so 
typically varies between 100 and 200 nm in InAs nanowires^^ 
and in the latter case E^o - 40 yueV only. 



Appendix F: Standard Rashba SOI and Rashba coefficient 

Both Ge and Si are inversion symmetric, and thus coupling 
of Dresselhaus type is absent. However, this does not ex- 
clude the conventional Rashba term (RSOI), Eq. (|7]). Here we 
briefly outline its derivation; details are described in Ref.lSOi 
As in Sec. lVIl we assume a constant electric field E^ along 
the X axis, which, referring to holes, results in the dipole term 
i/ed = —eExX as a perturbation added to the potential energy. 
Accordingly, H^d is added to the multiband Hamiltonian (en- 
velope function approximation), where it appears only on the 



diagonal, while off-diagonal parts provide the k ■ p coupling. 
Finally, a Schrieffer- Wolff transformation of the multiband 
Hamiltonian, with focus on the valence band Fg, yields the 
Rashba term 



Hso = aEJkyJ^ - k-Jy), 
^El 



(Fl) 
(F2) 



in third order of perturbation theory, where a is the Rashba 
coeflicient and additional, negligible terms have been omitted. 
In Eq. ( IF2b . £0 is the band gap (direct, k - Q) between con- 
duction (Fg) and valence (Fg) band, and P is the corresponding 
momentum matrix element between the i-like Fg and the p- 
like Fg,FiJ states^ For Ge, explicit values are £0 = 0.90 eV 
and P -9.1 eVA)^ which yields a x -0.4 nm^e. 

We can project Eq. dFlb onto the low-energy subspace 
{g+,g_,e+,e_) by calculating the 16 matrix elements. The ef- 
fective Hamiltonian for RSOI takes on the form 



//r = //™ = aE,S T^CT, + aE^k, ■ 



(F3) 



where 5 - {g+\kyj~\e+). This Hamiltonian has two effects: 
first, it features a constant coupling between the g and e states, 
and second, it provides a term which is linear in k^ and mixes 
the spin blocks. The latter is absent at k~ - 0, so that only 
the constant term aE^ST^cr- contributes for small k{, this is 
of the same form as the direct Rashba SOI //dr = eExUr^cr- 
(DRSOI) resulting from dipolar coupling. Finally, we note 
that 



aExS nm^ 



(F4) 



for Ge, so that the DRSOI dominates RSOI by one to two 
orders of magnitude in typical Ge/Si nanowires of 5-10 nm 
core radius. 



Appendix G: Coupling to magnetic fields 

In Eqs. (|9]) and ( fTOb . we show the effect of external mag- 
netic fields on the low-energy sector for fields applied along 
(z) and perpendicular (x) to the nanowire, respectively. Below, 
the explicit values for Z, and X, are listed. 



Zi = 0.75, 
Z2 = -0.81, 
Z3 = 23% R, 



X2 
X3 



= 2.72, 
= 0.17, 
= 8.04 /?, 



(Gl) 



using the parameters yi 
Ge^ 



13.35,7, = 5.11, and /c = 3.41 for 
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